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ABSTRACT 

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and 7-ray 
wavelengths of the BL Lacertae object AO 0235+164. We analyze the timing of multi-waveband 
variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) 
images at A = 7 mm with ^ 0.15 milliarcsecond resolution. The association of the events at different 
wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo 
simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum 
in the 7 mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering 
of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the 
quasi-stationary "core" and in the superluminal knot, both parsecs downstream of the supermassive 
black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated 
along the line of sight by light-travel time delays, that passes through a standing recollimation shock 
in the core and propagates down the jet to create the superluminal knot. The multi- wavelength light 
curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter 
time-scales. This, as well as the variability of the polarization and the dual location of the outburst, 
agrees with the expectations of a multi-zone emission model in which turbulence plays a major role 
in modulating the synchrotron and inverse Compton fiuxes. 

Subject headings: BL Lacertae objects: individual (AO 0235+164) — galaxies: active — galaxies: jets 
— gamma rays: general — polarization — radio continuum: galaxies 
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1. INTRODUCTION 

Our understanding of the processes leading to the gen- 
eration of 7-ray emission from blazars, the most extreme 
active galactic nuclei, depends on where those 7-rays 
originate. T his is currently the subjec t of considerable 
deba t e (e.g.. iMarscher fc JorstadI [20101 : iTavecchio et all 
[MntlAgudo et al.ll201in ~ Two main locations of the site 
of 7-ray emission in blazars have been proposed. The 
first is close (< 0.1—1 pc) to the supermassive black hole 
(BH) , which can readily explain short time-scales of vari- 
ability of a few hours reported in some 7-ray observations 
of bkzars (lAckermann et all 120101 : iFoschini et al.l 120101 
120111 : ITavecchio et al.ll2010f ). Although short time-scales 
of variability limit the size of the emission region and 
does not necessarily imply short distances to the BH, 
this scenario has the advantage that optical-UV photons 
from the broad emission-line region are available for scat- 
tering to 7-ray energies by highly relativistic electrons in 
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the jet. 

However, locating the 7-ray emission region so close 
to the BH contradicts the increasingly large number of 
time coincidences o f radio-millimeter and 7-ray events 
(IJorstad et alj 120101: IMarscher et"all 120101: lAgudo et alj 
1201 Ih . This problem is overcome if the 7-rays are 
emitted from a region much farther (>> Ipc) from 
the BH, beyond the "core" where the jet starts to be 
visible at millimeter wavelengths with very long base- 
line interferomctry (VLBI). Supporting this second sce- 
nario, Agudo ct al] (|2011[ ) unambiguously locate the 
region of 7-ray flares > 14 pc from the BH in the 
jet of OJ287 through correlation of millimeter-wave 
with 7-ray light curves and direct ultrahigh-resolution 
7 mm imagi ng with the V L BA. Simila r results are 
obtain ed by IMarscher et alJ (|2010D and IJorstad et alJ 
(I20T0I) for PKS 151 -089 and 3C 454.3, respectively. 
IMarscher fc JorstadI (|2010( ) have recently proposed a 
model that reconciles both 7-ray flare emitting regions 
located at >> 1 pc from the BH and intraday 7-ray vari- 
ability through a model involving turbulent plasma flow- 
ing through standing shocks in the core with a volume 
filling factor that decreases with energy and therefore 
frequency. 

This paper investigates the location and properties 
of a radio to 7-ray outburst in the BL Lacertae ob- 
ject AO 0235-hl64 (0235-^164 hereafter, z = 0.94). 
Historically, the blazar has exhibited extreme variabil- 
ity (by over an order of magnitude on time-scales 
< 1 year) across all spe ctral ranges, incl uding X-ray 
(e.g.. iRaiteri et all |2009() and 7-ray (e.g., lAbdo et"al] 
l2010al) ^ The most prominent radi o outbursts are ac- 
companied by op tical counterparts ( Ledden et al.l 119761 : 
iBalonek fc Dentlll980l : IRaiteri et al.l 12008'). although no 
correlation with variability at other wavebands has been 
reported previously. 

On sub-milliarcsecond scales, 0235-1-164 is extremely 
compact (5, 0.5 mas) at millimeter wavelengths (e.g., 
iPiner et al.l 120061) . and shows extreme superluminal ap- 
parent speeds (^^^'^ = 46.5±8.0c. ,Jorstad et al. 2001)^^ 
This speed, and th e variability Doppler factor i^ar = 
24.0 of 0235-^164 (jHovatta et all I2OO90 . sets a max- 
imum jet viewing angle d = arctan[2/3™p'^/(/3™p'^^ -f 
^lar - 1)] ^ 2?4. Although the large ^"^^^ value 
migh t lead to underestimation of d f or nq n-cylindric 
jets ()Gopal-Krishna. Wiita. fc Dhurdd 12006( 1 . the un- 
usual compactness of 0235-1-164 and its extreme flaring 
activity i s consistent with low Q values as estimated above 
(see also lHovatta et a"I]|2009() . Superluminal knots have 
been observed at position angles rangi ng from 5° to —55° 
(jJorstad et al.ll2001l : IPiner et al.ll2006( ). consistent with a 
broad (ttapp ~ 60°) projected jet opening-angle. Both 
Q and Ofapp constrain the maximum intrinsic jet half- 
opening-angle aint/2 = (ttapp sin^?)/2 < 1?25. 



2. OBSERVATIONS 

Our photo-polarimetric monitoring observations of 
0235-^164 (Figs. [HSl) include (1) 7 mm images with 
the VLBA from the Boston University monthly blazar- 

■^^ We adopt the standard ACDM cosmology, with Hq=1\ km 
s~^ Mpc~^, £7jv/ = 0.27, and Ha. = 0.73, so that 1 mas corresponds 
to a projected Unear distance of 7.9 pc. 



monitoring program , (2) 3 mm observations with the 
IRAM 30m Telescope, and (3) optical measurements 
from the following telescopes: Calar Alto (2.2m Tele- 
scope, observations under the MAPCAT^^ program). 
Steward Observatory (2.3 and 1.54m Telescopes^^), Low- 
ell Observatory (1.83m Perkins Telescope), San Pedro 
Martir Observatory (0.84m Telescope), Crimean Astro- 
physical Observatory (0.7m Telescope), and St. Peters- 
burg State University (0.4m Telescope). Our total flux 
light curves (Fig. [2]) include data from the Fermi-LAT 
7-ray (0.1-200 GeV) and ^m/t-XRT X-ray (2.4-10 keV) 
observatories, available from the archives of these mis- 
sions, and RXTE at 2.4-10 keV. Optical i?-band fluxes 
come from the Tuorla Blazar Monitoring Program^'', 
the Yale University SMARTS program^^, and Maria 
Mitchell and Abastumani Observatories. Longer wave- 
length light-curves were acquired from the Submillimeter 
Array (SMA) at 850 ^im and 1 mm, the IRAM 30m Tele- 
scope at 1 mm, the Metsahovi 14m Telescope at 8 mm, 
and both the Owens Valley Radio Observatory (OVRO) 
40m Telescope Fermi Blazar Monitoring Program^^ and 
University of Michigan Radio Astronomy Observatory'^'^ 
(UMRAO) 26m Telescope at 2 cm. 

We followed data reduction procedures d e scribe d 
m previous studies: VLBA: 'Jorstad et al.' (20051); 
optica l polarimetric data: [jorstad ct al. (2010); IRAM 
data: lAgudo et all (E^e', 2010,): SM A: iGurweh et all 
(12001: Me tsah ovi: [Rrasranta eFall (I1998D: OVRO: 
IRichards et"an (I2010h : UMRAO: l Aller et al.l (fl985l): 



Swift: IJorstad et al.( (|Mol ): RXTE: IMarscher et all 
(I2010D: and fe rmi-LAT: IMarscher etHI (|2010( ): 
lAgudo et al] poTlh . For both Swift and RXTE data, 
only the 2.4-10 keV band is included in the X-ray light 
curve. In this energy band, where the intrinsic spectral 
curvature a nd photoelect r ic ab sorption at z = 0.524 
discussed by IRaiteri et al.l ()2006l and references therein) 
are not important, a single absorbed power law (with 
Galactic absorption corresponding to a neutral hydrogen 
column density of 0.9 x 10^^ cm^) fit the data adequately. 
To process the 0.1-200 GeV Fermi-LAT data, we used 
2-day binning and a single power- law spectral fit wit h 
photon index held fixed at 2.14 ()Abdo et al.l l2010b[) . 
This produces essentially the sam e light curv e as th e 
broken-power-law spectral model of lAbdo et al.l ()2010c| ) . 

3. ANALYSIS OF OBSERVATIONAL RESULTS 

3.1. Major Millimeter Flare in 2008 Related to a New 
Superluminal Knot 

Our 7 mm VLBA maps of 0235-fl64 (Fig. [T]) reveal 
compact total intensity structure that is well fitted by 
one or two circular Gaussian components at most of the 
observing epochs, where the core generally corresponds 
to the brighter one. Between mid-2008 and mid-2009, 
the images also contain a second emission region (Qs, the 
brightest jet feature detected thus far in this object) that 
propagates at a superluminal apparent speed < /3app >= 

|http : //web . bu . edu/bleizars/VLBApro j ect . html] 
■^^ http : //www. iaa. es/~iagudo/research/MAPCAT 
http : // j ames . as . arizona . edu/~psmith/Ferini 
http://users.utu.fi/kani/lm 
http : //www . astro . yale . edu/smarts/glast 
http : //www . astro . caltech.edu/ovroblazars/ 
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Fig. 1. — Sequence of 7 mm VLBA images of 0235+164 convolved with a FWHM = 0.15 mas circular Gaussian beam. Images in our 
program before 2008-06-12 and after 2009-07-27, containing only a single emission feature (i.e., the core), are not displayed. Contour 
levels represent total intensity (levels in factors of 2 from 0.4 to 51.2% plus 90.0% of peak= 4.93 Jy/beam), color scale indicates polarized 
intensity, and superimposed sticks show the orientation of x- 



(12.6 ± 1.2) c. The position angle of QsG h 150°,- 
180°] contrasts with previously reported VLBI position 
angles by up to ~ 180°, which is consistent with a sudden 
change of the jet ejection angle amplified by projection 
effects. 

If the motion of Qs is ballistic, it was coincident with 
the core in 2008.30 ± 0.08, near the start of an extreme 
millimeter-wave outburst (labeled 08mm in Fig- HI that 
peaked on 2008 October 10 with a flux density ~ 6.5 Jy 
at 3 mm. Figure [2] shows that radio and millimeter- wave 
outbursts in 2008 (08rad and 08mm) contain contributions 
from both the core and Qs, whose fluxes reached maxi- 
mum on 2008 October 20 and November 16, respectively. 
Their contemporaneous co-evolution suggests that the 
disturbance responsible for the ejection of Qs extended 
from the location of the core to Qs in the frame of the ob- 
serve r, which could have resulted from light-t ravel delays 
(e.g., iGomez et all 119971 : lAgudo et aLll2001h . Qs is the 
brighter 7-mm superluminal knot ever seen in 0235+164, 
and flares 08rad and 08mm are the only outbursts that oc- 
curred after the ejection of Qs. The rarity of such events 
strongly implies that they are physically related. 

The jet half-opening-angle of 0235+164 (aint/2 < 
1°25) and the average FWHM of the core measured 
from our 31 VLBA observing epochs in [2007,2010] 
((FWHMcorc) = (0.054+0.018) mas), constrain the 7 mm 
core to be at dcoic = 1.8(FWIIMcorc)/tanaint > 12 pc 
from the vertex of the jet cone. 

3.2. Contemporaneous Flares from "f-ray to Radio 
Wavelengths 

Figure [2] reveals that the 08rad and 08mm flares were 
accompanied by sharp optical, X~ray, and 7-ray counter- 
parts (08opt, 08x, and 08-y flares, respectively). Our for- 
mal light -curve correla t ion an alysis (Fig. |4]) -performed 
following lAgudo et al.l ()2011[ )- confirms the association 
of 7-ray variability with that at 2 cm, 8 mm, 1 mm, and 
optical wavelengths at > 99.7% confidence. The flux 
evolution of the VLBI core is also correlated with the 



7-ray light curve at > 99.7% confidence. Moreover, 
the evolution of the degree of optical linear polarization 
(Popt) and X-ray light curve are also correlated with the 
optical _R-band, 1 mm, and 2 cm light curves at > 99.7% 
confidence (Fig. [5]), further indicating that the extreme 
flaring activity revealed by our light curves is physically 
related at all wavebands from radio to 7-rays. 

There is, however, no common pattern to the discrete 
correlation function (DCF) at all spectral ranges. This 
implies that, although there is correlation on long time- 
scales (years), on short time-scales (< 2 months) the vari- 
ability pattern does not correspond as closely. This is the 
result of the intrinsic variability pattern rather than the 
irregular time sampling at some spectral ranges. 

The sharp systematic peaks in the DCFs involving the 
i?-band light curve aids in the identification of relative 
time delays across wavebands, as measured with regard 
to the first sharp DCFr,a2 peak. In this way, we find 
that, relative to the i?-band maximum, the peaks of the 
A2 = 2 cm, 8 mm, 1 mm, and 7 mm core flares in 2008 are 
delayed by ~ 60, ~ 45, ~ 40, and ~ 40 days, respectively, 
whereas the X-ray delay is ~ 25 days. Only the 7-ray 
variations lead those at i?-band, by ^ 10 days according 
to the peak in the DCF seen in Figure U) 

3.3. Correlated Variability of Linear Polarization 

Figure [3] reveals extremely high, variable optical po- 
larization, popt > 30%, during the sharp 08opt optical 
peaks (see also IHagen- Thorn et all [2OO81 . for the 2006- 
2007 outburst). Whereas the integrated millimeter- wave 
degree of linear polarization (pmm) and that of the 7 mm 
core remain at moderate levels ^5%, the polarization of 
Qs (Pmm.Qs) pcaks at the high value of ~ 16 % close to 
the time of the second sharp optical sub-flare. The coin- 
cidence of this sharp maximum of Pmm.Qs in the brightest 
superluminal feature ever detected in 0235+164 with the 
(1) high optical flux and polarization, (2) flares across 
the other spectral regimes, and (3) flare in the 7-mm 
VLBI core, implies that the ejection and propagation of 
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Fig. 2. — Left: Light curves of 0235-1-164 from 7-ray to millimeter wavelengths. X-ray photon index evolution from the Swift-XKT data 
is also plotted. Vertical dotted lines mark the three most prominent OSopt optical peaks. The yellow area represents the time of ejection 
of feature Qs within its uncertainty. RJD = Julian Date — 2400000.0. Right: Same as left panel for RJDe [54500, 55000]. 



Qs in 0235-1-164's jet is pliysically tied to tiie total flux 
and polarization variations from radio to 7-rays. 

On long time-scales (years), the linear polarization an- 
gle at both optical (xopt) and millimeter (xmm and Xmm ) 
wavelengths varies wildly, without a preferred orienta- 
tion or systematic common trend. However, during flare 
OSopt, Xopt maintains a stable orientation at (100 ± 20)°, 
whereas xSm is roughly perpendicular to this (~ 0°), 
as expected for a plane-perpendicular shock wave prop- 
agating to the south towards Qs. Owing to the large 
peak value of Qs, Pmm,Qs 16 %, one cannot explain the 
orthogonal optical-millimeter polarizations by opacity ef- 
fects. Instead, we propose that the optical polarization 
mainly arises in a conical shock associated with the 7- 
mm core, while the millimeter-wave polarization results 
from a propagating shock front associated with Qs. We 
surmise that the moving shock also emits polarized op- 
tical radiation, since the optical polarization drops pre- 
cipitously when the orthogonal polarization of Qs peaks 
(FigMright). 



3.4. Low Probability of Chance Coincidences 

The relationship among the 7-ray, optical, and radio- 
millimeter flares is supported by probability arguments. 
If the flares occur randomly, the probability of observ- 
ing, at any time, a 7-ray outburst like the one reported 
here (i.e., with flux > 10~^ phot cm~^ s~^ and dura- 
tion ~ 70 days) is = 0.08. For optical and radio- 
millimeter wavelengths, this probability is Popt = 0.04 
and Pmm = 0.15, respectively. Thus, if the flares at 
different wavelengths were random and independent of 
each other, the probability of observing a 7-ray, op- 
tical, and radio-millimeter flare at any given time is 
P7,opt,mm = 5 X lO^"*. This couutcrs the null hypothesis 
of random coincidence at 99.95% confidence. This con- 
fidence keeps > 99.6 % even when unrealistically consid- 
ering twice as long as observed 7-ray, optical, and radio- 
millimeter flares, thus pointing out the robustness of this 
result. 



4. DISCUSSION AND CONCLUSIONS 
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The coincidence of the ejection and propagation of 
Qs — by far the brightest non-core feature ever reported 
in 0235+164 — with the prominent 7-ray to radio out- 
bursts and the extremely high values of Popt and Pmm.Qs 
provides convincing evidence that all these events are 
physically connected. This is supported by probabil- 
ity arguments and by our formal DCF analysis, which 
unambiguously confirms the relation of the 7-ray out- 
burst in late 2008 with those in the optical, millimeter- 
wave (including the 7-mm VLBI core) and radio regimes. 
Comprehensive studies of other blazars pla ces the core 
parsecs downstrearn of the BH Q orstad e t al.l 120101 : 
iMarscher et"al1 [20081 [2OIOI : I Agudo et al..,2011D . 

The general correlation of 7-ray and optical light 
curves (Figs. |4] and [5]) contrasts with the poor de- 
tailed correspondence during the main outburst (Fig. [5]). 
This is difficult to reconcile with the external Comp- 
ton mechanism, in which 7-rays arise from inverse 
Compton scattering of photo ns originating outside the 
jet by electrons in the jet (jBegelman fc Sikoral UMA 
iDermer fc Schhckeiseii 119931 ) ■ According to this model, 
rapid fluctuations of the 7-ray and optical flux are 
caused by changes in the number of electrons with suf- 
ficient energy to radiate at these frequencies. How- 
ever, such changes should affect the fiux at both wave- 
bands in a similar manner if the seed photon field varies 
smoothly with distance down the jet, as expected if 
these photon s are from the broad e mission-line region 
jSikora ct al. 1994)_ o r a d ust torus (jBlazeiowski et al] 
[2000; Ma lmros c c t aIll201lD . This is contrary to our ob- 
servations unless the magnetic field varies in a much dif- 
ferent manner than does the energy density of relativistic 
electrons. Instead, the nature of the variations suggests 
that the seed photon field also varies quite rapidly. 

A naturally variable seed photon source is synchrotron 
radiation from the jet. The maximum apparent 7-ray 
luminosity of 0235-1-16 4 is similar to that of the optical 
synchrotron radiation (jAbdo et al.l l2010al) . Thus, syn- 
chrotron self-Compton (SSC) radiation can potentially 
explain the 7-ray emission without leading to higher- 
order scattering that flattens the 7-ray spectrum. We 
can reconcile the short time-scales of variability, < 1 day, 
with a distance from the BH of parsecs, by appealing 
to both the very narrow opening angle of the jet (see 
above) and the possibility of structure across the jet. 
The former limits the width of the emission region to 
~ 0.8 pc, while the latter allows individual zones of emis- 
sion to be no larger than ~ 0.8/ \^ pc, where N is the 
number of zo nes, which presumab l y rep resent turbulent 
cells. In fact. IMarscher fc Jorstad (pOlOl ) used 0235+164 
as a blazar whose spectral energy distribution, polariza- 
tion, and variability properties can be well represented 
by such a model. The synchrotron flux varies in response 
to changes in the energy density of the flow, maximum 
electron energy, and magnetic fleld direction in each cell, 
while the SSC flux varies from the first two of these plus 
the time-delayed synchrotron photon density from the 
other cells. These variations follow similar trends, lead- 
ing to a general correlation but with detailed fiuctuations 
that can differ significantly. 

We identify the 7-mm core as the first re-collimation 
shock near the end of the jet's acceleration and colli- 
mation zone (ACZ I Jorstad et all 120071: IMarscher et"all 



I2OO8I I2OIOD . Superluminal feature Qs is consistent with 
a moving shock oriented transverse to the jet axis, given 
the extremely high Pmm,Qs7 with xSm parallel to the di- 
rection of propagation of Qs. The flux evolution of the 
core appears closely tied to that of Qs, and its light curve 
is correlated at high confidence with those at 7-ray, op- 
tical, and millimeter wavelengths. This suggests that Qs 
is the head of an extended disturbance, perhaps contain- 
ing a front-back structure stretch ed by light-travel delays 
in the observer's frame (see, e.g.. lAloy et al.l l2003'). The 
presence of this dual structure is also consistent with the 
perpendicular Xopt and xSm if corresponds to the 
front region emitting at miUimeter wavelengths, while 
the optical fiare arises mainly from the interaction of the 
slower back region with the standing conical shock in the 
core. 

Under this scenario, the radio/millimeter- wave and 
optical (and perhaps X-ray) synchrotron flares start 
when the front region crosses the conical shock at the 
core, where the jet is at least partially optically thin. 
This interaction accelerates electrons, increasing the syn- 
chrotron emissivity, and produces SSC 7-ray emission 
from the up-scattering of IR-optical photons. There is 
a time delay of the latter, since the seed photons must 
travel across part of the jet befor e arriving at the scatter- 
ing site fsee lSokolov et al.ll2004D . When the back region 
of the moving perturbation encounters the core, their in- 
teraction again produces efficient particle acceleration, 
which is seen as a sudden optical and radio/millimeter 
synchrotron emission enhancement. The subsequent op- 
tical variability is produced by the passage of the re- 
maining shocked turbulent plasma in the back structure 
through the core. During the different optical sub-flares, 
the integrated radio/millimeter synchrotron flux keeps 
rising. This radio/millimeter outburst is more prolonged 
owing to the longer synchrotron cooling-time of electrons 
radiating at these wavelengths and the lower speed of 
the back structure. Indeed, Qs does not reach its max- 
imum radio/millimeter- wave flux until traveling a pro- 
jected distance of ^ 0.13 mas from the core. 

When the entire front-back structure passes across the 
core, the synchrotron emission declines rapidly at opti- 
cal (and, if relevant. X-ray) frequencies, as does the 7-ray 
SSC emission. The decay of the radio/millimeter- wave 
emission is more gradual (see above) . The prominent X- 
ray flare on RJD^54760-54770 is consistent with being 
simultaneous with the peak of a millimeter and/or optical 
synchrotron flare and could arise from either scattering 
of millimeter photons or synchrotron radiation from the 
highest-energy electrons. The former is consistent with 
the evolution of the X-ray spectral index, which gradu- 
ally steepens to values consistent with the optically thin 
millimeter-IR spectral index (^ 1) at the time of the 
X-ray peak. 
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